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ABSTRACT 

The measurement of the brightness temperature fluctuations of neutral hydrogen 21 
cm lines from the Epoch of Reionisation (EoR) is expected to be a powerful tool for 
revealing the reionisation process. We study the 21 cm cross-correlation with Cosmic 
Microwave Background (CMB) temperature anisotropics, focusing on the effect of the 
patchy reionisation. We calculate, up to second order, the angular power spectrum 
of the cross-correlation between 21 cm fluctuations and the CMB kinetic Sunyaev- 
Zel'dovich effect (kSZ) from the EoR, using an analytical reionisation model. We 
show that the kSZ and the 21 cm fluctuations are anti-correlated on the scale corre- 
sponding to the typical size of an ionised bubble at the observed redshift of the 21 cm 
fluctuations. The amplitude of the angular power spectrum of the cross-correlation 
depends on the fluctuations of the ionised fraction. Especially, in a highly inhomoge- 
neous reionisation model, the amplitude reaches the order of 100 ^K^ at £ ~ 3000. 
We also show that second order terms may help in distinguishing between reionisation 
histories. 

Key words: cosmology: theory - cosmic microwave background - large-scale structure 
of the universe 
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1 INTRODUCTION 

The Epoch of Reionisation (EoR) is an essential milestone in the formation and evolution of cosmic structure. The first 
luminous objects produced in collapsed dark matter halos in the early universe {z ~ 20) started to reionise the inter galactic 
medium (IGM) which was neutral after recombination. Currently we have only a few observations for the EoR. The first one 
is the L y-g absorption m easurement towards high redshift QSOs which probes the fraction of neutral hy drogen along the line 
of sight ( Fan et al. 20061 ) . and the second one is the large-scale CMB polarisation ( Komatsu et al. 2O10l ). These observations 
indicate that the IGM was fully ionised by red shift z ^ 6. The recen t HST observations fou nd large samples of Lyman break 
galaxies (LBGs) at high redshifts, 7 < z < 10 i Bouwens et al. 20ld l. Bouwens et al. ( 2O10l ) have studied reionisation with a 
galaxy model based on these data. Their results suggested that, in addition to such high redshift LBGs, other reionisation 
sources, for example, faint galaxies and population III stars, are required to match the optical depth of the WMAP seven-year 
data. 

While, current observational data for the EoR are insufficient to study the details of the EoR. In recent years, several 
observations of signals from the EoR have been suggeste d to obtain furth er information about th e EoR, for example fluctuations 



of the neutral hydrogen 21 cm line (iMadau et al 



due to th e kinetic Sunyaev-Zel 'dovich (kSZ; 
review see 



1997 , for a review see Furlanetto et al 



20061). s mall-s c ale CMB anisot ropics 



Sunvaev fc Zel'dovichl ll98Q : lOstriker fc Vishniad 1 19861: IVishniad IiQStI. for 



. .Aghanim et al. 20081 ). and Ly-a damping of high redshift QSOs and gamma ray bursts I Miralda-Escude 19981 : 

Barkana fc LoeblboOll ). While the latter can provide us with information about the end of the EoR, the former two are expected 
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to probe the IGM during the EoR. The LOFAE0, MWA|1 and SKA|1 are being instaUed or designed for the measurement of 
21 cm hne fluctuations, while telescopes such as AClQ, SPlB and OLIMPO jMasi et al-lboOSf ) will be used to detect and 
measure the kSZ signal. 

Although both auto-correlations of 21 cm lines and CMB anisotropics during the EoR are good probes of the EoR, 
the cross-correlation between 21 cm fluctuations and the CMB anisotropics created during the EoR is also expected to be 
useful to study the history of the EoR. The cross-correlation has a potential to provide additional information other than 
their respective auto-correlations. Besides, the cross-correlation decreases the statistic errors caused by the foreground and 
the systematic effects, as compared to their auto-correlation. There are several analyt i cal or num erical works about the 
cross- correlation between CMB and 21 cm fluctuations during the EoR. Alvarez et al. ( 20061 ) and Adshead fc Furlanettol 
1 20081 ) computed the expected signal on large scales {I ~ 100) by analytically calculating the cross-corre lation between 21 cm 
fluctuations and the CMB Doppler anisotropics in the linear regime of the cosmological perturbations. Tashiro et al. ( 2O10l ) 
studied the detectability of these signals by LOFAR, MWA and SKA. On small scales {I > 1000), bec ause the domi nant 
contributions of CMB anisotropics come from the kSZ effect due to the patchiness of the ionised medium, ICooravi (|2004l ) has 
partially studied the cross-correlation with kSZ anisotropics and the second order 21 cm fluctu ations in a simple re ionisation 
model. He has also investigated the higher order cross-correlation by calculating the bispectrum. ISlosar et al.l (|2007l ) have also 
done the study of the 21 cm cross-correlation with the CMB SZ effect which is caused by hot electrons in the first supernovae 
remnants during the EoR. Since reionisation is a complex physical process, numerical simulations pla y an impor t ant ro le in 
the studies of the 21 c m cross-correlation with CMB temperature anisotropics. Numerical works by iJelic et al. I (|2010| ) and 



Salvaterra et al.l (120051 ). focus especially on the small-scale cross-correl ation due to the patc hy r eionisation. Additiona lly, the 



21 cm cross-correlation with CMB polarisation has been calculated bv I Tashiro et al. I (|2008l ) and iDvorkin et alj I (|2009l ). 

In this paper, we study the cro ss-correlation between kSZ anisotropics and the second order 21 cm fluctuations during 
the EoR analytically. ICooravi l|2004r ) has studied this cross-correlation in the simple analytical reionisation model where the 
fluctuations of the ionisation fraction are linearly related to the density fluctuations. He concluded that the cross-correlation 
cannot appear due to the geometric cancellation occurring between the velocity and the density fluctuations. However, the 
kSZ effect depends strongly on the evolution of the ionisation bubbles, and numerical studies of the c ross-correlation betwee n 
21 cm and kSZ anisotropics also shows that patchy reionisat ion generates si g nals on small scales IjSalvaterra et al.ll2005l ). 
Therefore, we revisit this issue with the analytical model of McQuinn et al. ( 20051 ) which produces a reionisation history 
similar to that found in recent numerical simulations. 

The outline of our paper is the following. In Sec. II, we give the analytical form of the second order cross-correlation 
between kSZ anisotropics and 2 1 cm fluctuations. In Sec. Ill, we give a short description of the analytical reionisation model 
based on iMcQuinn et al.l (|2005l ). In Sec. IV, we show the angular power spectrum of the second order cross-correlation and 
we discuss the detectability in the case of the SKA sensitivity. Section V is devoted to the conclusions. Throughout the paper, 
we use the concordance cosmological parameters for a flat cosmological model, i.e. h = 0.73 {Ho — h x 100 km/s/Mpc), 
To = 2.725K, fib = 0.05, = 0.27 and erg = 0.9. 



2 THE SECOND ORDER CROSS-CORRELATION 

In this section, we calculate the angular power spectrum of the cross-correlation between 21 cm fluctuations and kSZ 
anisotropics during the EoR at the second order in the fluctuations. For simplicity, we assume that both fluctuation flelds are 
isotropic statistically. Under this assumption, the angular power spectrum of the cross-correlation C'l is given by 

/ »kSZ 21 1 _ c-D r-D ^ f-,\ 

where a^f^i af^m^ multipole components of the CMB temperature anisotropics and 21 cm fluctuations during 

the EoR. 

2.1 kSZ CMB anisotropies 

During the EoR, secondary CMB temperature anisotropies are caused by the kinetic SZ effect. Their expression is 

/'70 
dr) g{'n)n ■ v{rj,n), (2) 



^ http://www.lofar.org 

^ http://www.mwatelescope.org/ 

^ http://www.skatelescope.org 

■* http:/ /www. physics. princeton.edu/act/ 

^ http://pole.uchicago.edu/ 
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where v is the baryon velocity field, g(rf) is the visibility function at the conformal time ry, and the present value of the 
conformal time is 770. The visibility function is given by (7(77) = fe~^ where r is the optical depth of Thomson scattering from 
r] to today and f = arXin^ with or the cross section of Thomson scattering, Xi the ionised fraction, and tih the neutral 
hydrogen density (we ignore the ionisation of helium) . 

We can decompose Xi and rzn into the background and fluctuation values, 

nH=nii(l + 5), Xi^Xi{l + 5x), (3) 

where the symbols with a bar represent the background values. In Eq. ((3|, since we assume that the hydrogen density follows 
the dark matter density on scales much bigger than the baryonic Jeans length, 5 is the total matter density fluctuation fleld. 
Substituting Eq. (|3]) into Eq. (O , we obtain 



rksz(n) = -Tcmb j dr^ g{r^)n- v{ri,n) {1 + S{r^,n) + 5^{rj,n) + 5{ri,n)S^{rj,'h)). (4) 
We focus on the second order part in Eq. ([4]), which can be written in terms of the Fourier components of the fluctuations 

as 

^TkszH = -»r„„b Jdril ^ I ^^i}'^ ^ Hv, fc - k'mn, k') + 5^{^, k')) exp[»(r,o - r,)(n ■ fc)], (5) 

where we use the relation r = [rjo — rf)h, and we relate the velocity to 5 by the continuity equation in the cosmological linear 
perturbation theory 

fc • 

f ^i^^i-n^k), (6) 

where the dot represents the derivative with respect to ry. 

Our flnal aim is to obtain the angular power spectrum of the cross-correlation. Therefore, we consider the spherical 
harmonic decomposition of Eq. ([Sjl, a^rrf = / dfi5T]^sz{'fi)YP . The spherical harmonic coefficients of the kSZ are given by 

j3i,. r ^'ii 



£ m £ m' 
l'" "m"' m"" 



X AS"r"''""" {v)5{-n, k^mn, fc2) + 5.,{V, fc2))^^^^ j."(fc2r)FA" {ki)Y^" (k), (7) 
where we replaced fe and fc' by fci = fe — fc' and fc2 = fc', and 



tmm'm"m"'m"" ■/ , m+m' -m" +m"" 647r .£' + £" / 3(2i?' + 1) ^t' t'" 1 r — mm' —m" m" " m -I \ I o\ 

Au'iiii'"! —^i Y 4tt{2£"" + ]^) '^-'"'-'^"'"'""'-"000 Jcmbffl,??). (a) 
Here, Ciif^*,„ are the Clebsch-Gordan coefficients and MJ"]™^™^*"* are the integrals of quadruple spherical harmonics, 
^.T™"' - j dhYr^-{h)Yr^-^{n)Yr/-{n)Yr^^{n) 

- C-n™! / (2^2 + 1)(2£3 + 1)(2£4 + 1) ^i^iil' ^l2t' It ^i2l'li (n\ 

- \ ^1 167r2(2£l + l) ^mamim'^OOO "-^mam' -mj *-"000 ■ V' ) 



2.2 21 cm fluctuations 



The brightness temperature of the 21 cm line from a redshift z is given as in lMadau et al.l l| 19971 ) by 



T21 



■{Ts - Tcmb){z), 



(10) 



where Tcmb is the CMB temperature and Tg is the spin temperature given by the ratio of the number density of hydrogen in 
the excited state to that of hydrogen in the ground state. The optical depth for the 21 cm line absorption T21 is 

3c^hAioXunu 



T2l{z) = 



(11) 



16kiyl^TsH{z)' 

where ^10 is the Einstein A-coefficient, 1/21 is the frequency corresponding to the 21 cm wavelength and xh is the fraction 
of neutral hydrogen, which is written as a function of the ionised fraction Xi = 1 — xh. Note that we drop the redshift 
space distort ion by the peculiar velo city fluctuations of neutral hydrogen in Eq. pip , although this effect enhances the 21 cm 
fluctuations (jBharadwai fc Alill2004l ). 

Combining Eq. ^ with Eqs. (|10p and (lllf) . we can obtain the observed 21 cm fluctuations at the observed frequency i^. 
The second order fluctuations which we here focus on is given by 



5T2i{fi,u) 



drj 



d^k d^k 



(27r)3 (27r)» 



o-W^2i(r?, rj{zobs))To(z(ri))5{rj, k - k')5H{'n, k') exp[i(T7o -rj){n- fc)]. 



(12) 
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where VK2i(r;, is the spectral response function of the observation experiment, normahsed as J drjWiiijj , ri{z)) = 1 and 

centred at rj{z), the redshift 2obs is related to the frequency v 'as = V2\/{1 + ^obs), 5h = {xh — xh)/xh and To is a 
normalisation temperature factor given by 



23 



0.02 



0.15 



l + z 



-,1/2 



mK. 



(13) 



The spin temperature is determined by three couplings with CMB, IGM gas and Ly-a photons. In the EoR, L y-a photons 
emitted from ionising sources couple the spin temperature with the IGM gas temperature ( Ciardi fc Madaull2003l ). Meanwhile, 
since the IGM gas is heated up quickly by Ly-a and X-ray photons from stars and QSOs, the IGM gas temperature is much 
higher than the CMB temperature during reionisation. Therefore, we can assume Ts 3> Tcmb during the EoR in Eq. (|13p . 
Taking the harmonic decomposition, we obtain the spherical harmonic coefficients of the 21 cm fluctuations, 

d^ki f d^k2 



y' m' f" m" 



(27r)S 



{r,)S{r,,ki)5„{r,,k2)ji'ik^r)ji„(k2r)Y,T *(fei)F,7? -(fea) 



where 



B" 



iv) 



(14) 



(15) 



Here A/™^J^|™^ is the integral of triple spherical harmonics, 



= (-1)' 



(2£2 + l)(2^3 + l)^^,^3f 



47r(2£i + 1) 



mi'^OOO 



(16) 



where mi + m2 



ma. 



2.3 The cross-correlation 

The second order cross-correlation is given by substituting Eqs. (O and ([14]) into Eq. ([l}. We obtain 



EE E I'^i^^-m 



d*k2 



d^k[ 



)3 J (27r)3 J (2tt)» J (27r)3 



{5*{ri, kr)5:{v, k2)5{n' , k[){S{v', k'2) + , ^2))) A^^T^p '^'^ 
X Jr, (fcir)jr/ (fer) ^^i^j,. (fc^^r^y.f * (fei)F,f * (fe 



m — m-, — m-i 



(17) 



where r' — rjo ~ rj' and we use Sx ~ —Sh- 

Under the assumption that all fluctuation flelds are Gaussian, the Wick theorem breaks the ensemble average in Eq. (|17|l 
into components with (55), (S^Sx) and (SSx)- For the simplification of Eq. (|17|l . we assume that W2i{z) = 5{z — 2oba). This 
is a good approximation because, compared to the observed frequency, the spectral resolution is narrow (for example, the 
spectral resolution in the LOFAR case is less than 1 MHz while the observed frequency is about 150 MHz for Zohs ~ 10). We 
can simplify further by using the approximation for the integration of spherical Bessel functions with ^ 3> 1, 

_ 71 F{£/r) 



dr' / dk k'^F{k)ji{kr)ji{kr') 



drf-^Fik) 



(18) 



Finally, we can rewrite the cross-correlation as 



kSZ-21 _ ^ (2^2 + l)(2^i -I- 1) ^^l^i2t^2 To{T]a],s)Tculh Gjrjobs) _ 



'^1^2 



2712(2^+1) 



girjobs) / dkji^{krobs) 



dr 



£2 



Psx rjobs, — + Pxx Vobs, — PssiVobs, k) + [ Pss { Vobs, — + Psx { V, — ] ] Psx{rjobs, k) 



where the power spectra Pss, Pxx and Psx are defined as {5{r],ki)5{rj,k2)) = (27i)''5(fci — k2)Pss{'il,k\), {5x{r),k\)5x{'rj,k2)) = 
(27r)^(5(fci - k2)Pxx{v,ki), and {5{f], ki)5x{ri, k2)) = (27i)^5(fci - k2)Psx{v, ki). In Eq. (US]), G is the growth factor of the dark 
matter density fluctuations which is 5{k,r]) = G{ri)5{k) with the present density fluctuations 5{k). Now, the epoch we are 
interested in is matter dominated, so that we can assume Goc 1/(1-1- z) in terms of the redshift 2. 

In order to calculate the cross-correlation, the power spectra Pxx and Psx which are determined by the reionisation model 
are essential. We discuss the analytical reionisation model in the following section. 
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3 REIONISATION MODEL 



For an analytical reionisation model, we adopt the approach of lFurlanetto et alj l|2004bf ) and lMcQuinn et alj (|2005l ). lonisation 
bubbles start to evo lve from high density galaxy regions into the voids, as shown in recent numerical simulations (e.g. 
Trac fc Gnediiilboogl . and references therein). Therefore, the mass of ionised gas mion is associated with the mass of a collapsed 
object rrigai by the Ansatz, mion = C™-gai where C, is an ionizing efficiency. The condition for the full ionisation of a region of 
mass m is that the region contains sufficie nt sources to self- ion ise, i.e. /coii > C~^i where /coii is the fraction of collapsed halos 
above the critical mass for collapse, iLacev fc Colelll993l ). 



This criterion gives the barrier (the density threshold) 5x for "self-ionisation" which depends on m. iFurlanetto et al 



(2004a|) found a reasonable approximation of the barrier in the linear form of the variance of the density fluctuations, 
a^{m,z), as B{m,z) = Bo + Bia^{m,z) where a{m,z) is obtained by smoothing the density field at the scale m. Here, 
Bo^Sc- V2K{(:)amin{z) and Bi 
for the collapsed ionisation source 
For the linear barrier B(m, z) 



- d5x/da^\„2^Q where CTmin(z) is the mass dispersion at the minimum mass and redshift z 
the bubble mass function is written as I Sheth 19981 ) 



dnim] 
dm 



dm - 



2 p 



dlogcr 



dlog m 



Bo 



exp 



B^{m, 



2a^{m) 



dm , 



(20) 



where p is the mean mass density of the Universe. 

The smallest bubble mass is given by i^mmin. Therefore, we can obtain the mean ionised fraction (volume averaged) Xi as 



/ 



dn{m) 
V(m] — ^ — -dm - 
dm 



-2SoSi 



erfc 



Bo 



+ ^erfc 



Bo + Bicr? 



(21) 



where o"^ = (j{C,m, z) and V{m) is the comoving volume of a bubble with mass m. 

In the case of a linear barrier, the linear bias of a source of mass m is given by (jMcQuinn et aLlboOSl ) 



6(m) 



1 + 



B{m)/a^{m) - 1/Bo 



D{z) 



Therefore, the mean bias of the bubble fo(m) is obtained from 



dm b{m)V{m' 



dn{m) 
dm 



(22) 



(23) 



In this reionisation model, the free parameters for the model are C, and mmin. Here we ta ke two parameter sets which 
are motivated from numerical simulations: "stars" model and "QSOs" model ( Jelic et al. 2O10l ). In both models, the ionised 
fraction reaches Xi = 0.5 ai z — 11, in order to agree with the WMAP results. 

In the "stars" model, we assume that stars are responsible for reionisation. We take a low efficiency ^ ~ 40 which is 
reasonable for normal star formation and assume that the minimum mass corresponds to a virial temperature of 10^ K, above 
which cooling by atomic hydrogen becomes efficient. 

In the "QSOs" model, we assume that the reionisation history is faster and the bubble size is larger compared to those 
in the "stars" model. Therefore, we set high a virial temperature (5 x 10'' K) and a high efficiency « 200. The candidates 
for the ionisation sources are massive stars and QSOs. 

We show the evolution of the ionised fraction for each model in Fig. [T] From Eq. H20p . we can obtain the bubble size 
distribution Vdn/dR as a function of the comoving size i? of a bubble under the assumption that the bubbles are spherical. 
We plot the results in Fig. (2] 



3.1 The two-point correlation function (,xx{r) 



In order to obtain the power spectra P^x and Pgx in Eq. p9|) . we need to compute the correlation function (,xx{r) = 
{xi{xi)xi{x2)) — Xi and S,sxir) ~ (5i xi)xi(x2)) where the p oints Xi and X2 axe separated by r = \xi — X2\. Here we utilize 
the analytical correlation functions of iMcOuinn et al.l ( 20051 ). 

As in the case of the density correlation function in the halo formalism, the correlation function of the ionised fraction 
^xx{r) receives two contributions. One is a one bubble term Pi which is the two-point correlation for the case where two 
points which are separated by r are ionised by the one and same ionisation source, the other is a two bubble term P2 which 
corresponds to the case where two points are ionised by two separate sources. As shown in Fig. [2] the typical size of an 
ionisation bubble becomes larger than 5 Mpc when the ionised fraction reache s one half In such regi me, where the ionisation 
bubbles become large. Pi is largely dominant and P2 can be ignored. Thus, McQuinn et al. 1 20051 ) divide the reionisation 
process into two phases: the early phase and the late phase. In the early phase, both Pi and P2 are important, while in the 
late phase. Pi is dominant and P2 can be ignored. The criterion for these phases is set as Xj > 0.5 in order t o be in agreement 
with results from the hybrid approach of analytic modeling and numerical simulations of lZahn et al. I l|2005l ). They define the 
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Figure 1. Evolution of the mean ionised fraction. The soUd and dotted hnes represent Xi in the "stars" and "QSOs" models, respectively. 
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Figure 2. Ionised bubble comoving size distribution. The dotted, solid and dashed lines represent the distributions in the "stars" model 
a,t z = 10, z = 11 and z = 12, respectively. The ionised fractions are Xi = 0.8 at z = 10, Xi = 0.5 at 2 = 11 and Xi = 0.3 at z = 12. We 
also plot the distribution in the "QSOs" model at 2: = 11 as the thin solid line. The left side of each line ends at R(C"*min) where C™min 
is the minimum mass of the ionised region. 



correlation function ^xx{r) by 



where 
Pi(r) = 

P2{r) = 



(1-x, 
Pi(r)4 



A(r) 
P2(r) 



when Xi > 0.5, 
otherwise, 



dm '•^"'('^^ 
am 

dn{m\ 

ami ; 

am 



r), 



f'dm,^^^ r 

dm J 



d r2[l + ^(ri — r2\mi,m2) 



(24) 

(25) 
(26) 



Here, (,{r\mi,m2) is the excess probabihty to have a bubble of mass mi at the distance r from a bubble of mass 7112. For the 
simplicity of the calculation, it is assumed that ^{r\mi, 7712) can be written in terms of the correlation function of the matter 
density ^ss as 5(r|mi, 7712) — 6fiA(max(r-, Ri + R2)) where i?i(mi) and ^2(1^2) are the bubble radii. 

In order to calculate the volume in Eqs. (|25p and (|26p analytically, all ionisation bubbles are assumed spherical. Therefore, 
Vo{m,r) is the volume within a sphere o f mass m that can encompass two points separated by a distance r. For the volume 



integration in Eq. (|26p . iMcQuinn et al.l (|2005f ) adopt the overlapping conditions; (1) mi cannot ionize r2, and m2 cannot 



ionize ri ; (2) the center of m2 cannot lie inside mi , but any other part of m2 can touch mi , and vice versa. 
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3.2 The two-point cross-correlation function £,5x('i') 

As in the case of 5ix(r), the two-point cross-correlation £,sx{'r) has two contributions, Pin and Pout. The contribution Pin 
corresponds to the case of both points being contained within the same ionised bubble. Following iMcQuinn et al.l (|2005l ). it 
is written as 

■mu dnh{mh\m) 



Pn(r) = I dm^^^Vo{m,r) J dnm 



p drUh 



dm^^^4^M'm;r)[l + B{m,z)], (27) 
dm 



where the last line in Eq. (I27p is obtained by using the fact that the inner integral is the mean over-density of the bubble 
1 + Sb and 5s is B{m, z) at linear order. 



The contribution Pout corresponds to the case when one point is outside the ionised bubble of the other point. iMcQuinn et al, 

1 2005h give Pout in terms of the mean bias for halos bh, 

Poutir) = J dm^^^^Vo{m,r) + J dm^^^ J d'r4bHHss{r~ri,)]. (28) 

where dnh{mh\m) /dnih is the conditional mass function. In Eq. (|28p . the integration range of is over all bubbles which ionise 
the point but not the other point separated by r from rb. For simplicity, £^ss is evaluated at the separation max[P(m),r]. 

As the reionisation proceeds and the typical size of an ionised bubble becomes large, the term Pout becomes unimportant 
as compared to Pi„. Therefore, the computation of is divided into two phases again, 

Pin -Pi when 0.5, 

1 Pin + Pout — Xi otherwise, 

where we assume that Pn is dominant in large Xi (xi > 0.5) and we subtract Pi given by Eq. (|27p from Pn which is the 
correlation between Xi and p/p. 



4 RESULTS AND DISCUSSION 

We calculate the angular power spectrum of the cross-correlation described in Eq. p9|) in the two models, "stars" and 
"QSOs". First, we show the results in the "stars" model in Fig. [S] In this model, the mean ionised fraction is 0.3, 0.5 and 
0.8 at z = 12, z = 11 and z = 10, respectively. The signal of the cro ss-correlation between kSZ and 21 cm fluctuations 



exhibits an anti-correlation on small scales {£ > 1000). As mentioned in ICooravl (|2004l ). there is a geometric cancellation in 
the cross-correlation. This cancellation is responsible for a suppression of the amplitude of the cross-correlation. However, 
the cross-correlation has a distinctive oscillatory shape. Especially, we found that the peak position of the anti-correlation 
represents the typical size of an ionised bubble at each redshift. For example, at z — 11, the typical size of an ionised bubble 
is almost 6 Mpc, as shown in Fig. [2l and the anti-correlation at Zobs = 11 is maximal at the corresponding multipole £ ~ 4000. 
As the Universe evolves, the typical scale of an ionised bubble becomes larger. The peak position of the anti-correlation shifts 
accordingly toward smaller values. 

The evolution of the cross-correlation amplitude depends on the evolution of through the power spectra of Pxx and P^s 
which evolve rapidly during the EoR. Since the amplitudes of Pxx and PxS increase as the redshift decreases, the amplitude 
of the cross-correlation also becomes larger at low redshifts. However, after the average ionisation rate reaches Xi ~ 0.9, the 
signal of the 21 cm fluctuations becomes weak and the cross-correlation amplitude also starts to decrease. 

In Fig. [31 we also plot the firs t order cross-correlat ion between 21 cm and CMB Doppler anisotropics calculated by using 



the same expression as Eq. (15) of lAlvarez erZI (|200d '). The sign of the first order cross-correlation depends on the evolution 



of 5x- As long as Sx is small, the ionisation process is homogeneous, and the cross-correlation is negative. On the other hand, 
in the case of a highly inhomogeneous reionisation, the sign of the cross-correlation is positive. In our reionisation model, 
the first order cross-correlation at the early phase of reionisation is negative at £ < 1000 (see the top and middle panels in 
Fig. O. We found the amplitude of the first order cross-correlation at Zobs ~ H is 300 pK^ at the peak position, £ ~ 100, 
and decreases rapidly towards zero at large multipoles. As we can see in Fig. [31 the second order kSZ-21 cm cross-correlation 
dominates the first order cross-correlation at multipoles larger than £ — 1000. However, as the ionisation process proceeds, the 
ionisation fraction is highly inhomogeneous and Sx is evolved well. As a result, the first order cross-correlation has a positive 
sign and a high amplitude as shown in the bottom panel of Fig. [3l The first order cross-correlation becomes comparable to 
the second order kSZ-21 cm cross-correlation even at £ ~ 1000, while the kSZ cross-correlation still dominate the first order 
cross-correlation and has negative correlation at multipoles higher than £ = 1000. 

Next we show the dependence of the angular cross-correlation power spectrum on the ionisation model in Fig. |4l In the 
"QSOs" model, the ionisation history is rapid and the typical size of ionised bubbles is large. The amplitude of Pxx and PxS 
in the "QSOs" model is larger than in the "stars" model. As a result, in the "QSOs" model, the signal of the cross-correlation 
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Figure 3. Angular power spectra of the second order cross-correlation in the "stars" model. From top to bottom panels, we plot the 
angular power spectra at Zq^s = 12, Zobs = H and Zq^s = 10, respectively. The mean ionised fraction is x = 0.3 at z = 12, x = 0.5 
at z = 11 and x = 0.8 at z = 10. For reference, we show the first order cross-correlation between the CMB temperature and 21 cm 
fluctuations as the dotted line in each panel. 



is large and the peak position of the anti-correlation appears on small multipoles, as expected. We can therefore conclude 
that the cross-correlation between kSZ and 21 cm fluctuations at the second order is sensitive to the average size of an ionised 
bubble. The first order cross-correlation also has a higher amplitude than in the "stars" model because the amplitude depends 
on the evolution rate of the backgroun d ionis ation fraction. However, the inhomogeneous contribution coming from the term 
with PxS in Eq. (15) in Alvarez et al. ( 20061 ) is partially canceled by the homogeneous one from the term with Pss- As a 
result, in the highly inhomogeneous "QSOs" reionisation model, the cross-correlation between kSZ and second order 21 cm 
fluctuations reaches a significant amplitude, compared with the first order cross-correlation at small scales (I < 1000). 



4.1 Detectability 

In the previous section, we showed that the peak position of the anti-correlation is related to the typical bubble size at the 
observed redshift of 21 cm fluctuations. Here, our concern is the detectability of such negative peak in the kSZ-21 cm cross- 
correlation. In order to investigate the detectability, we calculate the signal-to- noise ratio (S/N). For simplicity, we assume 
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Figure 4. Dependence of the cross-correlation on tlie ionisation model. The solid and the dashed lines are the power spectrum for the 
"stars" model and the "QSOs" model, respectively (see text). We set ziobs = H where x = 0.5 in both models. For reference, we plot the 
first order cross-correlation in each model as the thin lines. 



that CMB, 21 cm fluctuations and instrumental noise are Gaussian. The total S/N can be calculated as 

^max |^21-CMB|2 

/sky ^ (2^ + 1) I „2i-CMB 12 I n^incyin ' 



]\l J Jsky ^ 1 |(^21-CMB|2 (^21(^CMB ' 

where /sky is the sky fraction common to the two cross-correlated signals, and Cf^^, Cf^ and (^^i-cmb ^^.^ ^j^^ angular 
power spectra of CMB, 21 cm fluctuation and the cross-correlation between 21 cm and CMB, respectively. In order to focus 
on the detectability of the signal from the typical bubble size, we set ^min = 500 and ^max ~ 5000. 

At the multipoles that we are interested in {£ > 1000), the dominant CMB signal is due to the thermal SZ effect 
i Zel'dovich fc Sunvaeviri969 l). 'lowever we can remove this contribution because of the frequency dependence of the SZ effect. 
Therefore with the assumption that the foreground can be completely removed from the CMB map, the main contribution 
to Cp^^ comes from the primordial CMB anisotropies C^"^' and the noise of the instrument Np^^ . We can write Cp^^ as 

^« = Cf exp(-^VSMB/2) + iVf ^^ (31) 
where we assume the beam proflle of CMB observation is Gaussian with the Full Width at Half Maximum of the beam 9cmb, 
and o"cMB = ^CMB / VSIn 2. The effect of the beam size is a da mping of the signal of the primordial CMB on smaller scales 
than the FWHM. The noise power spectrum Np'^^ is given by (jKnoxlllQQsl ) 

where (Jpix is the sensitivity in each pixel and Slpix is the solid angle per pixel; fipix = ^cmb- 

As for the 21 cm fluctuations, the noise signal from the instruments and foreground will dominate the intrinsic signal 
from the EoR. Assuming that the foreground can be removed to the level below the noise from instruments, we can write 



27r 



2 



tobsAi. v^/ry^ ^^^^ 

where we use the noise power spectrum of 21 cm observation estimated bv lZaldarriaga et al.l (|2004l ). In Eq. psp . Ai/ is the 
bandwidth, fobs is the total integration time, A/T is the sensitivity (an effective area divided by the system temperature) and 
D is the length of the baseline associated with the FWHM of the 21 cm observation 621 = X/D. 

In the calculation of the cross-correlation signal, we assume that the foregrounds and noise of 21 cm fluctuations and CMB 
anisotropy are not correlated. Therefore, the cross-correlation consists mainly of the first order Doppler-21 cm cross-correlation 
and the second order kSZ-21 cm one, 

^21-CMB|2 /|^21 — Dopncri2 . 1^21— kSZi2\ r n2 / 2 , 2 \ /nl /oa\ 

\Ct \ ={\Ci "^"^ I +\Ci j )exp[-£ (crcMB +cr2i)/2], (34) 



where a2i ~ S2i/V81n 2 and both signals are affected by the angular resolution of the observations. 

Our interest is the detectability of the cross-correlation signal from the patchy reionisation by Planck and SKA. Therefore, 
in the computation of Eq. H31|) . we adopt the typical value of Planck which are ^cmb — 5 arrcmin and cTpix = 5 X 10"*^. The 
goal sensitivity of SKA is currently designed as A/T = 5000 m'^K^^ at 200 MHz. The configuration area is 20 % of total 
collecting area for 1 km baseline, 50 % for 5 km baseline, 75 % for 150 km baseline. Because we are interested in the scales 
£ ~ 2000, we take D = 1 km and A/T — 1000 m^K~^. The sky fraction /sky corresponds to the one of SKA because we 
consider Planck as CMB observation, which is almost full-sky. We assume 200deg^ per field of view and 4 independent survey 
fields for SKA. Therefore the total sky fraction is /sky ~ 0.02. 
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Figure 5. The S/N ratio for the detection of the cross-correlation signal at z = 11 as a function of the observation time. The left panel 
is for the "stars" model and the right panel is for the "QSOs" model. In both panels, the solid and dotted lines represent S/N for SKA 
and for the observation v^ith a 5 times better sensitivity than that of SKA, respectively. We set /gjjy ^ 0.02 in all plots. 



We plot S/N as a function of tots in units of hours for "stars" and "QSOs" models at z = 11 in Fig. [5| In both panels in 
Fig. [5] S/N of SKA with Planck is represented by the sohd hues. Obviously, longer observation times make S/N larger. Then, 
since the cross-correlation amplitude in the "QSOs" model is higher than in the "stars" model, S/N in the former model is 
lager than in the later model. However both S/N are below the detection level. This difficulty of the detection is mainly due 
to the instrumental noise of the 21 cm observation. Although the primary CMB is one of the si gnificant sources of noise in the 
detec tion of the cross-correlation signal between CMB and 21 cm from EoR on large scales ( Jelic et al. 20ld : Tashiro et al 



2O10l ). the primary CMB suffers Silk damping on the scales we are interested in here and the noise of Planck is also kept 



below the sufficient level. 

In order to clarify the impact of the improvement in the sensitivity of 21 cm observation, we calculate S/N in the case 
of a 5 times better sensitivity than that of SKA and plot the result as the dotted line in Fig. [5| The improvement of the 
sensitivity of 21 cm observation brings large S/N. Especially, the S/N in the "QSOs" model can reach S/N ~ 5 in 500-hour 
observation. Finally, while we use the same sky fraction /sky ~ 0.02 in all calculations, larger sky fractions also make 5'/A'' 
higher. 



5 CONCLUSION 



We investigated the small scale cross-correlation between CMB anisotropics and the 21 cm fluctuations during the EoR 
in harmonic space. The CMB anisotropies at small scales are mainly caused by the kSZ effect which is the second order 
fluctuation effect generated by the peculiar velocity and the fluctuations of the visibility function. We therefore calculated the 
cross-correlation with the second order fluctuations of 21 cm fluctuations. 

The cross-correlation signal between kSZ an d 21 cm fluctuations is negat ive on smal l scale s. This anti-correlation on small 
scales was found in the numerical simulations of Salvaterra et al. (2005) and Jelic et al. ( 2O10h . We found that the position of 
the negative peak is at the angular scale corresponding to the typical size of an ionised bubble at the redshift probed by 21 cm 
fluctuation measurements. This angular scale shifts to larger scales as ionised bubbles evolve. The amplitude also increases 
with the reionisation process until the average ionisation fraction reaches Xi ~ 0.9. The amplitude of the cross-correlation 
strongly depends on the typical bubble size. The cross-correlation in the case of larger bubbles has a higher amplitude than 
in the case of smaller bubbles, even if in both cases the mean ionisation fractions are the same. Moreover, the amplitude of 
the cross-correlation from large ionised bubbles is comparable to that of the first order cross-correlation. Those characteristic 
features of the cross-correlation could be used to distinguish between different reionisation histories with future observations. 

We also estimated the detectability of the small-scale cross-correlation by the current design sensitivity of SKA. It is rather 
difficult to detect the cross-correlation signal even in the radical reionisation cases. However, if the sensitivity is improved by 
a factor of 5, the detection or non-detection of the cross-correlation signal will definitely provide information about the EoR. 



ACKNOWLEDGMENTS 

HT is supported by the Belgian Federal Office for Scientific, Technical and Cultural Affairs through the Interuniversity 
Attraction Pole P6/11. 



© 0000 RAS, MNRAS 000, 000-000 



Second order cross-correlation between kSZ and 21 cm fluctuations from the EoR 11 



REFERENCES 

Adshead P. J., Furlanetto S. R., 2008, MNRAS, 384, 291 

Aghanim N., Majumdar S., Silk J., 2008, Reports on Progress in Physics, 71, 066902 

Alvarez M. A., Komatsu E., Dore O., Shapiro P. R., 2006, Astrophys. J., 647, 840 

Barkana R., Loeb A., 2004, ApJ, 601, 64 

Bharadwaj S., Ah S. S., 2004, MNRAS, 352, 142 

Bouwnes R. J., et al., 2010, arXiv: 1006.4360 

Ciardi B., Madau P., 2003, Astrophys. J., 596, 1 

Cooray A., 2004, Phys. Rev. D, 70, 063509 

Dvorkin C, Hu W., Smith K. M., 2009, Phys. Rev. D, 79, 107302 

Fan X., et al., 2006, Astronomical Journal, 132, 117 

Furlanetto S. R., Oh S. P., Briggs F. H., 2006, Physics Report, 433, 181 

Furlanetto S. R., Zaldarriaga M., Hernquist L., 2004b, Astrophys. J., 613, 1 

Furlanetto S. R., Zaldarriaga M., Hernquist L., 2004a, Astrophys. J., 613, 16 

Jelic v., et al., 2010, MNRAS, 402, 2279 

Knox L., 1995, Phys.Rev.D, 52, 4307 

Komatsu E., et al., 2011, ApJS, 192, 18 

Lacey C, Cole S., 1993, MNRAS, 262, 627 

Madau P., Meiksin A., Rees M. J., 1997, Astrophys. J., 475, 429 

Masi S., et al., 2008, Memorie della Societa Astronomica Italiana, 79, 887 

McQuinn M., Furlanetto S. R., Hernquist L., Zahn O., Zaldarriaga M., 2005, Astrophys. J., 630, 643 

Miralda-Escude J., 1998, ApJ, 501, 15 

Ostriker J. P., Vishniac E. T., 1986, ApJL, 306, L51 

Salvaterra R., Ciardi B., Ferrara A., Baccigalupi C, 2005, MNRAS, 360, 1063 

Sheth R. K., 1998, MNRAS, 300, 1057 

Slosar A., Cooray A., Silk J. I., 2007, MNRAS, 377, 168 

Sunyaev R. A., Zel'dovich I. B., 1980, MNRAS, 190, 413 

Tashiro H., Aghanim N., Langer M., Douspis M., Zaroubi S., 2008, MNRAS, 389, 469 
Tashiro H., Aghanim N., Langer M., Douspis M., Zaroubi S., Jelic V., 2010, MNRAS, 402, 2617 
Trac H., Gnedin N. Y., 2009. larXiv:0906.4348l 
Vishniac E. T., 1987, ApJ, 322, 597 

Zahn O., Zaldarriaga M., Hernquist L., McQuinn M., 2005, Astrophys. J., 630, 657 
Zaldarriaga, M. and Furlanetto, S. R. and Hernquist, L., 2004, Astrophys. J., 608, 622 
Zel'dovich Y. B., Sunyaev R. A., 1969, Astrophys. Space Sci., 4, 301 



© 0000 RAS, MNRAS 000, 000-000 



